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We find that magnetic neutrino-electron scattering is unaffected by oscillations for vacuum mixing 
of Dirac neutrinos with only diagonal moments and for Majorana neutrinos with two flavors. For 
MSW mixing, these cases again obtain, though the effective moments can depend on the neutrino 
energy. Thus, e.g., the magnetic moments measured with Ve from a reactor and from the Sun 
could be different. With minimal assumptions, we find a new limit on using the 825-days 
SuperKamiokande solar neutrino data: < 1.5 x 10~^''^b at 90% CL, comparable to the existing 
reactor limit. 



In the minimally-extended Standard Model, neutrinos 
of mass nil, have tiny loop-induced magnetic moments 
/iy ~ 3 X lO^^^/iB (m-iy/1 cV), where //s is the Bohr 
magneton. In various extensions of the Standard Model, 
larger magnetic moments can occur without large neu- 
trino masses. In the presence of flavor mixing, the funda- 
mental magnetic moments are associated with the mass 
eigenstates (since either a boost or a magnetic moment 
can be used to reverse the helicity). In the mass eigen- 
state basis, Dirac neutrinos can have diagonal or off- 
diagonal (transition) moments, while Majorana neutrinos 
can only have transition moments 

In the current experiments, the effects of neutrino mag- 
netic moments can be searched for only in the recoil 
electron spectrum from neutrino-electron scattering . 
Below we consider the interplay between magnetic mo- 
ments and flavor mixing for this process. We show how 
magnetic moments can be defined for beams that are 
initially neutrino flavor eigenstates. In some important 
cases these moments do not oscillate, i.e., they do not 
depend on distance from the source. However, in the 
presence of MSW mixing, these defined flavor moments 
can differ from the vacuum case and can depend on the 
neutrino energy, though not on the distance. As an illus- 
tration, we derive a new limit on the magnetic moment 
from the SuperKamiokande (SK) solar neutrino data [||. 

There are two incoherent contributions to neutrino- 
electron scattering: weak scattering, which preserves the 
neutrino helicity, and magnetic scattering, which reverses 
it. Thus the differential cross section is given by 
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In Eq. (0), QL = sin^ 9w + 1/2 for v^, Ql = sin^ 9w - 1/2 
for and Vn and — src? 9w for all flavors (for an- 
tineutrinos, exchange ql and gji). The magnetic moment 
pLu is expressed in units of ^b- Magnetic scattering, the 
second term in Eq. (|l]), grows rapidly with decreasing 
electron recoil kinetic energy T. 



In principle, there can be weak-magnetic interference 
effects. There is a negligible effect due to the fact that a 
massive neutrino is not a helicity eigenstate Also, if 
the neutrinos have a transverse polarization, the electron 
azimuthal angle distribution can be affected ; we ignore 
this case, as the effects are presently unobservable. 

Vacuum Mixing: The effects of flavor mixing on 
the weak scattering are well-known. Whatever the com- 
position of the neutrino beam, the different flavors are in 
principle distinguishable and hence their cross sections 
combine incoherently, weighted by the probabilities for 
the neutrino to be of each given flavor. 

We want to explore how neutrino oscillations affect 
the magnetic scattering. The shape of the electron recoil 
spectrum in magnetic scattering is universal (the same for 
all mass eigenstates) . The only quantity that depends on 
the beam composition is the effective magnetic moment 
^y. Let us assume that we begin with a beam of electron 
neutrinos. Under the usual oscillation hypothesis such a 
beam propagates over the distance L from its source in 
vacuum according to 
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where Uek is an element of the unitary mixing matrix and 
k labels the mass eigenstates. 

Similarly to above, whatever the composition of the 
neutrino beam, the different mass eigenstates are in 
principle distinguishable in the magnetic scattering, and 
hence their cross sections combine incoherently, weighted 
by the squares of the amplitudes for the neutrino to be of 
each mass after the scattering. Then the combined cross 
section for magnetic scattering has the form of Eq. (|^) 
with magnetic moment squared fi^ given by 
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where the summations j, fc, k' are over the mass eigen- 
states, and the subscript e labels the initial flavor. We 



1 



have made the usual relativistic expansion and have de- 
fined the oscillation length Lkk' = 47ri?y/Am^j,, for k ^ 
k' (there is no L dependent phase for k — k'). The quan- 
tities /ijfe in Eq. (|^) are the fundamental constants (in 
units oi hb) that characterize the coupling of the neu- 
trino mass eigenstates to the electromagnetic field. The 
summation over j is outside the square because of the in- 
coherence of the cross sections for different final masses. 
The expression for /ig simplifies in some important cases. 

Let us assume first that the neutrinos are Dirac par- 
ticles (with n flavors) with only diagonal magnetic mo- 
ments {^jk = Hj6jk)', this is the scenario used by the 
Particle Data Group §|. Then 



(4) 



and there is no dependence on the distance L or neu- 
trino energy E,^. In this case one can characterize the 
magnetic scattering by the initial flavor index instead 
of the mass indices. (Hence we disagree with Ref. in 
which the magnetic scattering depends on the final flavor 
index, i.e., it oscillates). Measurements of all magnetic 
moments and mixing parameters would allow extraction 
of the "fundamental" moments fij . 

Next consider the case of Majorana neutrinos, and as- 
sume that only two mass eigenstates are relevant. Then 

M? = |Ml2p(|C/elP + |C/e2n = |Ml2p, (5) 

which is not only independent of the source distance and 
the neutrino energy, but also of the mixing angle. 

Under what circumstances does one have to worry 
about a dependence on distance and neutrino energy, 
and in particular, when can such terms be dominant? 
Clearly, at least one term of the type /ijk x /Ujfc' with 
k ^ k' must be nonvanishing and as large as or /Lt|j,, . 
In other words, in the 3x3 matrix iijk there should be 
at least two comparable entries on the same line (and in 
the same column). For the Dirac case this implies that at 
least one nondiagonal magnetic moment is as large as the 
diagonal ones. For the Majorana case it implies that two 
different nondiagonal magnetic moments are of a similar 
magnitude. Both of these cases seem unnatural. 

MSW Mixing: The above discussion must be mod- 
ifled for matter-enhanced oscillations (the MSW effect). 
First, the initial composition of the beam is governed 
not by the vacuum mixing angle 9y, but by the initial 
matter mixing angle dm, which depends on lS.vr? j^En 
and the electron density. If the initial density is well 
above the resonance density, as is true for the standard 
solutions to the solar neutrino problem, then d„i "^1:11 
to an excellent approximation. Then initially, = 
cos6',„|j/i) -I- sin6',„|i^2) — \v-i)- 

Second, although a nearly pure \v-2) is produced in the 
solar center, if the passage through the resonance is nona- 
diabatic, then the final beam can be a mixture of and 
\v-2). Most generally [fo|, the mass eigenstates evolve as 
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where |cip -I- |c2p = 1. The phases 0a and 06 (real func- 
tions that depend on integrals of the instantaneous mass 
basis eigenvalues) are irrelevant here, due to the non- 
interference of different mass eigenstates in the magnetic 
scattering. For the adiabatic case (e.g., the solar large- 
angle solution ||ll|), C2 = 0. For the nonadiabatic case 
(e.g., the solar small-angle solution jl^), and a narrow 
resonance region (which naturally obtains), the probabil- 
ity of hopping from one mass eigenstate to the other is 
Pho-p = |c2p, which depends on the neutrino energy but 
not the distance from the source, e.g., for an exponential 
density with density scale height rg. 
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Thus for two-flavor Dirac mixing with only diagonal 
moments, we obtain for the effective magnetic moment 
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Note that this is different from Eq. (Q), even in the adi- 
abatic case. However, for the two-flavor Majorana case, 
we again obtain — |/ii2p, as in Eq. In both cases, 
since the initial state is a pure |i'2), there are no interfer- 
ence terms that depend on distance. 

SK Data: The best direct limit on the neutrino mag- 
netic moment, 1.8 x 10"^° fi^ at 90% CL comes from 
studies of neutrino-electron scattering with reactor an- 
tineutrinos. (See Ref. and references therein for the 
astrophysical limits). As explained above, the meaning 
of the measured /i^ using solar neutrinos and reactor an- 
tineutrinos could in principle be different. Nevertheless, 
it is important to realize that a magnetic moment nu- 
merically equal to the current reactor limit would have a 
statistically significant effect on the solar neutrino data 
from SK. Since there is, as explained below, no evidence 
in the data for a nonvanishing magnetic moment, we de- 
rive, with a minimum of assumptions, a limit on what we 
caU 

If the expected weak scattering rate were known (as as- 
sumed in Ref. jl^), an observed excess in the total rate 
would indicate a nonzero magnetic moment. However, 
as the total weak rate is not known a priori, we instead 
look at the shape of the electron spectrum for the effects 
of a magnetic moment. The signature of a nonvanishing 
magnetic moment would be an enhancement, compared 
to the weak scattering alone, of the events at low recoil 
energies, with less enhancement at higher energies. That 
is not observed. Instead, as shown below, the electron 
spectra recorded by SK have, within the statistics, the 
shape that one expects from weak scattering alone (we 
show below that the deviations observed currently in the 
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FIG. 1. The scatter plot illustrating the ranges and fre- 
quencies with which the fitted a and fj.^ appear. The plot 
corresponds to Qre/ = 0.5 and /ire/ ~ 0, and uses the relative 
errors (504-days data) of Ref. M. 



highest energy bins are irrelevant for our purpose). How- 
ever, the total number of events is less than the standard 
solar model predicts, presumably due to neutrino oscilla- 
tions. We do not need to know the value, or the mixing 
mechanism behind it, of this overall reduction of the scat- 
tering rate. 

We assume only that the shape of the measured spec- 
trum is not due to a fortuitous cancellation between a 
magnetic moment effect rising at low energies and an 
oscillation effect rising at high energies. The Sudbury 
Neutrino Observatory will check the spectral shape and 
total flux of the component. 

The procedure we adopt uses the measured relative er- 
rors by SK and the fact that the measured shape agrees 
with expectations. We calculate {da/dT) by folding 
Eq. (|l|) with the neutrino spectrum from Ref. [p^ . 
For both weak and magnetic scattering, we include the 
SK energy resolution though it makes little differ- 
ence in the final results. We histogram the results in 0.5 
MeV bins in total electron energy, as in SK. Thus, as a 
function of the bin number i, we have constructed the ex- 
pected spectra nw{i) and UM^i) for weak and magnetic 
scattering, respectively. 

In order to determine the upper limit of |/ie°'P we must 
take into account the statistical fluctuations in the SK 
data. While the data points divided by the solar model 
expectation are consistent with an energy independent 
reduction factor a, the individual bins are distributed, 
presumably randomly, around that value. To take that 
into account we choose some reference values aref and 
fj,ref and create a set of simulated data, ns{i), which are 
Gaussian-distributed around the theoretical expectation 
arefnw{i) + fJ'lef''^M{'i) with the relative errors a{i) given 
by SK. 

We then minimize the x^, 
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FIG. 2. The upper panel shows the relative frequencies 
of different fitted /il values for the indicated values of /ire/ 
and for Oref = 0.5. The calculation uses the relative errors 
(504-days data) of Ref. [^. The lower panel allows one to 
find for any fitted /i^ the most probable value of the true /ij 
(full line) and the 10% and 90% CL interval of that quantity 
(dashed lines). 
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with respect to the fit parameters a and /i^. For fixed 
aref and /i^^j , we repeat this procedure many times and 
plot the frequencies with which given values of the fitted 
a and /x^ appear. An example of the scatter plot of the 
fit parameters is shown in Fig.|l|. One can see, naturally, 
that the most probable values of the fit are the refer- 
ence values aref and ^ref- Also, the two variables are 
strongly anticorrelated (correlation coefficient r ~ —0.9), 
i.e., larger a is accompanied by smaller iif,. Dividing the 
numerator and denominator in Eq. ( |lOf ) by nw{i)^ one 
sees that the depends only on = ns{i)/nw{i), i.e., 
precisely on the quantities pubhshed by SK 

By repeating the calculation for different values of /i^e/ 
and projecting on the /i^ axis, one gets the distributions 
shown in the upper panel of Fig. ^. These distributions 
are Gaussian, and their width is almost independent of 
the chosen value of /i,-e/. Based on them we obtain the 
lines in the lower panel of Fig. ^ signifying confidence 
levels at 10%, 50% (the mean), and 90%. For a given 
fitted /i^ obtained in an experiment, these allow one to 
determine the likely range of the true /i^. For example, 
if one found a fitted //^ ~ 4, then from Fig. ||, the most 
probable true value is /i^ ~ 4, with the upper limit being 
~ 8 and the lower limit being ~ 0. Similarly, for a fitted 
/j,^ ~ 0, the true //^ is ^ 3.9. This would be the largest 
true nl, , due to statistical fluctuations of the finite data, 
that could have given this fitted /i^ ~ 0. 

In this way we solve problems associated with the 
statistical fluctuations as well as with the constraint 
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TABLE I. Limits on the magnetic moment (at 90% CL, in 
units of 10~^^ hb) and statistical errors (la) on a for various 
data sets j^: 504 days (I); the same just up to 13.0 MeV, 
to exclude the bins with an apparent excess compared to the 
constant data/SSM ratio (II); 708 days; and 825 days. In all 
cases we ignore the bin from 14.0 - 20.0 MeV. 



case 


1 sol 1 2 


1 sol 1 

W \ 


Sa 


504 days, I 


< 3.9 


< 2.0 


0.025 


504 days, II 


< 4.2 


< 2.0 


0.027 


708 days 


< 2.5 


< 1.6 


0.018 


825 days 


< 2.3 


< 1.5 


0.017 



l^f'P > 0. Figure 2 was calculated with are/ — 0.5, 
as observed in SK, but doesn't change significantly for 
0.4 < Uref < 0.6. Note that the results summarized in 
Fig. ^ can be also obtained analytically, without gener- 
ating many simulated spectra. The conclusions, in par- 
ticular the lower panel of Fig. H, simply follow from the 
properties of the individual sums in Eq. (|l0|). 

Using the SK data the fitted a ~ 0.5; the exact 
value is irrelevant since we are testing only the spectral 
shape, and not the normalization. The fitted values of 
fil are slightly (but not significantly; see Fig. 2) nega- 
tive: ~ —5, —3, and —2 (in the units of Table I) for the 
504-, 708-, and 825-days data sets. The slight (but di- 
minishing with time) positive slope observed in the data 
cannot be caused by a magnetic moment (which causes 
an increase at low energies) , though it could be caused by 
oscillations. The most conservative conclusion is there- 
fore to say that the slope is not negative, i.e., that the 
fitted fi^ values are not positive. That is, we obtain the 
limit by using Fig. 2 (and its analogs) and an assumed 
fitted value of /i^ = 0. Thus the limits in Table I are 
slightly weaker than what is naively implied by the data, 
but are more robust. The sensitivity to improves 
with time only as t~^^^, but the addition of more low- 
energy bins (e.g., the two added since the 504-days data) 
gives a more dramatic improvement. 

The uncertainties Sa in Table I reflect the increase in 
the error in the parameter a when one allows magnetic 
scattering. Our procedure does not include the corre- 
lations between systematic errors in different bins and 
therefore will not reflect the full systematic uncertainty. 
Note that in the standard analysis Q one assumes /i^ = 
and hence the uncertainty Sa is reduced by the factor 
l/Vl ~ r'^ — 2.3. Similarly, if the value of a were accu- 
rately and independently known, and we fit for /i^ only, 
an identical improvement in the upper limit of |/ie°'P 
would result. 

Conclusions: In this paper, we present three new 
results. First, that while neutrino magnetic moments 
are most fundamentally defined for mass cigcnstatcs, in 
several cases of practical interest non-oscillating (i.e., in- 
dependent of distance) effective magnetic moments can 
be defined for the flavor eigenstates. For Dirac neutrinos 



with only diagonal moments, these results are Eq. (Q) for 
vacuum mixing and Eq. (^ for MSW mixing. For Majo- 
rana neutrinos with two flavors, the result is /Xg = |/xi2p, 
for either vacuum or MSW mixing. Second, that MSW 
mixing can change the definition of the effective magnetic 
moment (allowing a dependence on the neutrino energy) , 
so that the measured moments using D,, from a reactor 
and i/fi from the Sun could be different. Third, that the 
shape of the SK recoil electron spectrum can be used to 
place a limit on the neutrino magnetic moment (note that 
we do not invoke any mechanism for neutrino interaction 
with the solar magnetic field). In general, this is a new 
limit, independent of the limit from reactor studies (with 
the same meaning only for solar vacuum oscillations) . In 
any case, the limit obtained using the preliminary 825- 
days data, < 1.5 x 10~^°/zb, is comparable to the 
existing reactor limit of 1.8 x 10~^^iib fl^ - 
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